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Abstract 


A  formalism  is  proposed  for  treating  the  problem  of 
ionization  of  adsorbed  species  on  solid  surfaces.  The 
ionizing  agents  are  taken  to  be  impact  atoms  and  laser 
radiation  with  frequency  low  compared  to  the  inverse  of 
characteristic  collision  times.  The  physical  constraints 
of  short  collision  times  and  low  laser  frequency  then 
allow  one  to  treat  the  adatom-surface-plus-field  system 
under  the  quasi-static  approximation  (QSA)  and  the  impact- 
atom-adatom-surface  collision  dynamics  under  the  impulse 
approximation  (IMA) .  The  latter  leads  to  a  time-dependent 
ionization  cross-section  which  is  factorizable  into  the 
square  of  an  electron-atom  scattering  matrix  element  and 
a  spectral  function  describing  the  energy-momentum  distri¬ 
bution  of  electrons  in  the  adatom-surface-plus-field  system. 
The  formalism  focuses  on  the  spectral  function  which  is 
shown  to  be  derivable  from  a  single-particle  Green's  func¬ 
tion  exactly  calculable  for  the  present  problem. 


1 


I.  Introduction 

Photoemission  studies  of  adsorbed  species  on  solid  surfaces, 
both  theoretical1-3  and  experimental,4  6  have  produced  consider¬ 
able  information  on  the  nature  of  chemical  bonding  between 
adsorbed  atoms  or  molecules  and  solid  surfaces.  The  process  of 
collision-induced  ionization  of  surfaces,  of  potential  importance 
in  generating  the  same  information,  has,  however,  claimed  rela¬ 
tively  little  attention,  due  to  the  widespread  belief  that  it 
would  be  extremely  difficult  to  extract  useful  experimental 

information  from  such  studies.  In  this  area,  new  ground  has 

7 

been  broken  recently  by  Conrad  et  al.,  who  considered  theoret¬ 
ically  the  problem  of  surface  Penning  ionization  of  a  single  CO 
molecule  chemisorbed  on  a  Pd (111)  surface  by  metastable  He*-beams 
and  compared  their  results  with  existing  experimental  data.  On 
another  front,  gas-phase  studies  of  Penning  ionization  also  suggest 

that  laser  radiation  can  have  pronounced  effects  on  the  emitted- 
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electron  energy  spectrum  '  as  well  as  the  ionization  probability 
itself.1®  These  works  all  point  to  the  likely  fruitfulness  of 
carrying  out  studies  on  collision-induced  ionization  of  adsorbed 
species  on  solid  surfaces  in  the  presence  of  laser  radiation.  In 
this  paper,  we  will  propose  a  formalism  to  treat  this  problem. 

The  physical  situation  we  are  considering  may  entail  several 
competing  processes  in  addition  to  the  emission  of  electrons.  There 
may  be  desorption  and  migration  of  adatoms,  or  even  formation 
of  free  ions.  However,  the  degree  of  catastrophe  induced  on  the 
surface  by  the  atom  beam  and  the  laser  radiation  is  not  entirely 
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beyond  the  experimentalist's  control.  Laser  power  and  laser  fre¬ 
quency,  for  instance,  could  be  selected  to  minimize  desorption 
and  migration,  and  the  incident-atom  impact  energy  could  conceiv¬ 
ably  be  adjusted  such  that  free-ion  formation  does  not  compete 
significantly  with  pure  scattering.  Also  the  amount  of  internal 
excitation  carried  by  the  impact-atoms  may  be  tailored  to  prefer¬ 
entially  ionize  the  adsorbed  atoms  rather  than  those  of  the  bulk 
medium. 

The  main  advantage  of  the  laser  as  an  inducing  tool  in  the 
present  case  is  that  it  has  much  greater  versatility  here  than, 
say,  in  the  process  of  laser-induced  desorption,  because  no 
resonance  requirement  on  the  laser  frequency  need  be  imposed. 

We  will,  however,  assume  that  the  laser  frequency  is  much  less 
than  characteristic  band  structure  resonances  of  the  pure  metallic 
surface  so  that  photoemission  need  not  be  considered  as  a  competing 
process.  We  also  require  it  not  to  be  in  resonance  with  adatom 
vibration  or  phonon  coupling  modes  so  as  to  avoid  dealing  with  de¬ 
sorption  or  migration  of  adsorbed  species. 

Since  the  laser  frequency  is  considered  to  be  low,  most  of  the 
energy  required  for  ionization  will  have  to  be  supplied  by  the  in¬ 
ternal  excitation  energy  and  the  translational  kinetic  energy  of 
the  incident  atoms.  Hence  we  consider  projectile  atoms  in  an  excited 
state  with  large  impact  velocities,  leading  to  short  collision  times 
t  such  that  t  <<  u”1,  where  u>L  is  the  laser  frequency.  This  means 
that  the  laser  photons  will  only  have  a  relatively  short  time  in 
delivering  energy.  But  since  the  adsorbed  atom  (adatom)  states  can 
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be  considered  to  be  broadened  by  the  solid  surface  into  a  "near¬ 
continuum  of  states"*3  (with  large  uncertainties  in  orbital  energies)  , 
the  low  frequency  photons  may  still  be  effective  in  transferring 
energy  to  electronic  degrees  of  freedom  in  the  adatom-surface  system, 
facilitating  electron  "hops"  between  the  adatom  and  the  surface.*4 
This  would  not  be  the  case  if  sharp  resonance  electronic  states  were 
considered  (such  as  in  gas-phase  collisions) . 

The  conditions  of  short  collision  time  and  low  laser  frequency 
permits  the  use  of  the  quasi-static  approximation  (QSA)  where,  even 
though  the  total  Hamiltonian  is  time-dependent  (due  to  the  radiation 
interaction) ,  the  energy  of  the  system  is  considered  to  be  adiabatic- 
ally  conserved  within  the  duration  of  a  characteristic  collision 
time  (cf.  Eqs.(2.6)  and  {2.1)  below].  Also,  the  smallness  of  t  war¬ 
rants  the  use  of  the  impulse  approximation*5  (IMA)  in  the  treatment 
of  the  projectile  atom-adatom-surface  collision  dynamics.  Within 
this  approximation,  the  collision  between  the  projectile  atom  and 
the  adatom-surface  (AS)  system  is  assumed  to  be  mediated  by  a  single  elec¬ 
tron  possessing  a  characteristic  momentum  and  energy  distribution 
determined  by  virtue  of  its  being  part  of  the  adatom-surface  plus 
field  system  and  otherwise  considered  to  be  free.  Naturally  the  IMA 
will  be  more  suitable  when  applied  to  cases  where  the  adatom-surface 
system  to  be  ionized  has  loosely  bound  electrons.  The  momentum  and 
energy  distribution  will  be  most  conveniently  obtained  through  a 
Green's  function  formalism.  The  IMA  has  been  successfully  applied 

to  a  wide  variety  of  collision  processes,  such  as  fast  electron-atom 

16  17  18 

collisions,  A(p,2p)B  scattering  in  nuclear  reactions,  '  and  gas- 

19 

phase  collisional  ionization.  Recently,  we  have  also  applied  it  to 

20 

high-energy  positron  ionization  of  adsorbed  species,  a  process 
closely  related  to  that  discussed  in  the  present  work. 
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In  what  follows  we  will  construct  a  formalism  for  the  calcu¬ 
lation  of  the  differential  cross  section  of  the  ionization  process. 
This  formalism  is  based  largely  on  many-body  techniques  leading  to 
the  construction  of  an  adatom-surface-plus-field  Green's  function 
that  is  time-dependent.  The  implementation  of  the  main  approxima¬ 
tions,  the  QSA  and  the  IMA,  will  be  shown  in  the  course  of  the  de¬ 
velopment. 


II.  The  Ionization  Cross-Section 

We  consider  the  process  in  which  projectile  atoms  B  with 
momentum  p^  are  incident  on  an  adatom-surface  (AS)  system  which  is 
driven  by  monochromatic  laser  radiation  with  field  strength  repre¬ 
sented  classically  by  E(t)  =  EgCOsaj^t,  where  is  the  laser  fre¬ 
quency.  Assuming  that  this  collisional  process  leads  to  ionization 
with  emitted  electron  momentum  and  final  momentum  for  the  atoms  B 

equal  to  p  and  pf  respectively,  the  time-dependent  differential 

17 

ionization  cross-section  can  in  general  be  written  in  the  form 


d3pd3p 

d6a(t)  =  ~  £|T(t)  |25{e+e.+E(+)  (t)-(E  (t)+e.+A)}i - 1 

*(+>  f  0  1  v0  ( 2uh)  J 


(2.1) 


In  Eq. (2.1)  vn  is  the  incident  velocity  of  B  atoms; 


e  =  p  /2m, 


(2.2) 


ei  =  Pi/2lV 


(2.3) 


ef  =  Pf/2mB 


(2.4) 


are  the  kinetic  energies  of  the  emitted  electron  and  the  free  atoms 
B  before  and  after  the  collision  respectively,  with  m  =  mass  of 


f.  -,L  s  % 
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electron  and  =  mass  of  atom  B;  and 

A  =  Ei  -  Ef,  (2.5) 

the  internal  energy  transfer,  is  the  difference  between  the  initial 
and  final  internal  energies  of  the  atom  B.  The  presence  of  the  6- 
function  with  time-dependent  energy  variables  follows  from  the 
quasi-static  approximation  (QSA) ,  where  it  is  assumed  that  if  the 
collision  time  t  is  short  enough  compared  with  the  period  of  the 
driving  force,  the  ground  state  energy  Eg  and  excited  energies  E^ 
of  the  unperturbed  and  singly-ionized  AS+field  system  vary  adiabatic- 
ally  over  t  <  x;  also,  the  total  energy  of  the  system  is  adiabatically 
conserved.  The  QSA  then  implies  that 

“AS  (t)  I  ^0  (t)  >  ~  E0<t)W0(t)>,  (2.6) 

and 

HAS(t)  |^(+)  (t)>  a  E(+)  (t)  |f'(+)  (t)>,  (2.7) 

where  HAg (t)  is  the  total  time-dependent  Hamiltonian  of  the  AS  system 
plus  field;  and  |H'Q(t)>  and  [¥^(t)>  are  the  corresponding  adiabatic 
ground  state  and  excited  state  wave  functions,  unperturbed  and  singly- 
ionized,  respectively.  The  summation  in  Eq.(2.1)  is  over  all  the 
final  singly-ionized  AS  states  |'^+*(0)>,  where  it  is  assumed  that 
at  t=0,  the  laser  field  is  turned  on.  Within  the  impulse  approxima¬ 
tion  (IMA) ,  the  transition  matrix  element  T  can  be  written  as 


T(t)  -  <y<+>(t)|a(p1)|l'0(t)>TEfi5f>5jEi>5ii5ii 


(2.8) 


where  a(p^)  is  a  fermion  annihilation  field  operator  (removing  an 
electron  of  momentum  p^  from  the  AS  system  in  its  ground  state  at 

is  an  off-shell  electron-atom  scattering 


time  t)  and  T 


E£,P£,p?E^,p^,p^ 


matrix  element,  describing  the  collision  between  an  atom  B  in  the 


initial  state  | i>  with  kinetic  energy  p?/2mB  and  a  free  electron 

with  momentum  to  produce  an  atom  B  in  the  state  | f >  with  kinetic 

2 

energy  p^/2m^  and  a  free  electron  with  momentum  p.  The  IMA  to  T  is 

schematically  represented  in  Fig  1,  in  which  and  p^  represent  the 

energy  and  momentum  removed  from  the  AS+ field  system  at  time  t. 

The  IMA  also  implies  that  at  the  right  vertex  (blank  circle)  energy 


and  momentum  are  both  strictly  conserved.  Hence  we  have 

p? 


and 


el  +  211^  +  A  ~  2^  +  2m  ' 


P1  +  Pi  =  P  +  Pf‘ 


(2.9) 


(2.10) 


We  first  discuss  the  sum  in  Eq.(2.1).  This  equation,  together 
with  Eqs . (2.8)  and  (2.9),  implies  that  the  sum  can  be  written  as: 

S(p1,e1,t)  =  J  |<H>(+)  (t)  |a(px)  !H'0(t)>|26(e]+E(+)  (t)-EQ(t)) 


=  l  <Vn(t)  laf(p.  )&(e.+E{'n  (t)-En(t))  |'F'"r;  (t)> 

(+)  0  u 

x  (t)  |a(p1)  |¥0(t)> 


,(+) 


(+) 


*  ^<Y0(t)  |aT(p1)6(£1+HAS(t)-E0(t))U^+)  (t,0)  |  V  K^’  (0  )  > 

X  <'F(  +  )  (0)  |u^+)  +  (t,0)a(p1)  |V0  (t)> 
»  <y0 (t) |a+(p1)6(E1+HAS(t)-E0(t))a(p1) |Y0(t)>  ,  (2.11) 


(  +  ) 


where  in  the  second  equality  use  has  been  made  of  Eq.(2.7)  (the  QSA) , 
and  the  time  evolution  operator,  U^+'>  (t,0)  ,  for  the  Schrddinger  picture 
^  (t)  has  been  introduced  such  that 


wave  function  V 


«r  ■  • 


Ug  (t,o)  I'Fq  (o)>  =  ! (t) >  • 


(2.12) 


The  quantity  S(p^,e^,t)  is  referred  to  as  the  time-dependent 
spectral  function  and  can  be  interpreted  as  the  probability  at 
time  t  of  finding  an  electron  with  momentum  p^  in  the  AS+ field 
system  and  the  energy  E^tJ-e.^  in  the  residual  system  after  an 
electron  has  been  removed.  The  differential  ionization  cross- 
section  can  then  be  written  as: 

3  ’ 

db0(t)  =  S(p.,e  ,t)  jT  D  D  •  E  t>  D  '  V - 1  *  (2*13) 

n  ix  /Pf /Pyt-^Pj^P^  vQ  ^2.  j 

The  time-dependent  spectral  function  can  be  shown  be  related 
to  the  Fourier  Transform  of  the  advanced  single-partic  ^reen's 
function  as  follows.  We  start  with  the  definition  of  the  advanced 


single-particle  Green's  function; 


G£  (j?i ;  t ,  t ' )  =  i<H'0|aH(p1,t)aH(p1,t')  |H'o>0(t-t’) 


(2.14) 


where 


V  2  IV0)>' 


(2.15) 


and  0(t-t')  is  the  Heaviside  step  function.  aH(p^,t)  is  the  time- 
dependent  Heisenberg  picture  fermion  annihilation  field  operator 


given  as: 


aH  (^1 '  fc>  =  ' 


(2.16) 


with  Ug ( t , 0 )  satisfying  the  relation 
|  Y0  (t)  >  =  Ug(t,0)  |<F0>. 

We  can  then  write  the  advanced  Green's  function  as 


(2.17) 
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GA(pl;t»t')  = 


i<'F0|Ug(tf0)a+{p1)Us(t,0)Ug(t'  ,0)a(p1)Us(t*  ,0)  |  VQ>e  (t-f  ) 


=  i<4'0|Ug(t,0)a+(p1)Us(t,t,)a(p1)Us(t+t1,t)Us(t,0)  |4'0>6(-t1) 
=  i<4'0|Ug(t,0)a+(p1)Ug(t+t1,t)a(p1)Us(t+t1,t)Us(t,0)  |  4<0>  9  (- 

=  i<¥0(t) |a+(p1)Ug(t+t1,t)a(p1)Ug(t+t1,t) |¥0(t)>9 (-tj.) . 


(2.181 


In  the  second  equality  above  we  have  made  the  change  of  time  variables 


4  "  f  -  t 


(2.19) 


and  use  of  the  property  for  the  time  evolution  operator  Ug  that 


*  vw- 


(2,20) 


Under  the  QSA  we  can  write 

Us(t1+t,t)  *  exp{-  i  HAS(t)t1} 


(2.21) 


if  1 1^ |  <  t  and  it  follows  that 


ga  (pi ;  t ,  t+t  i )  =  i<H'0(t)ja  (p-, )  e 


^(HAg(t)-Eo(t))ti 


GC-^Jatpj^)  |'?(J(t)>, 


t1\  <  t.  (2.22) 


Introducing  the  integral  representation  for  the  step  function 
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GA(pl'a,/t)  = 


<Vt} ! a*  (Pl> 


a(px)  |4'0(t)>, 


(2.25) 


From  Eq.(2.11),  the  relation  between  the  spectral  function  and  the 
restricted  advanced  Green's  function  is  then 

S(P1,e1/t)  =  limfG^p^o^-in/t)  -  Gft  ,a)1+in  ,  t)  ] 

n+o 


5S  rm  VPi-V1"'*1- 


(2.26) 


where 


and  the  representation  for  the  delta  function 

X  _  1 _  lim , _ 1  _ 1  v 

0K  ’  2ttT  n+o  vx^Tn  x+Tiv 


(2.27) 


(2.28) 


has  been  used.  We  note  that  in  arriving  at  Eq.(2.25)  the  restricted 
Fourier  Transform  is  taken  as  though  Ug(t^+t,t)  assumes  for  all  t-^  the 
functional  form  dictated  by  1 ^ I  <  x,  that  is,  the  Hamiltonian  enterina 

-h  A 

into  the  computation  of  GA(p^;t,t’)  is  just  HAS (t)  [cf.  Eq. (2.21)1. 

The  restricted  transform  is,  of  course,  distinct  from  the  true  trans¬ 
form 


GA(pl'“,t)  =  c^tl  e^^^^G^  (Pi ;  t ,  t  ’ )  ,  t^=t'-t. 


(2.29) 


The  computation  of  GA(p^,t,t')  will  be  discussed  in  the  next  section, 

^  « 

and  that  of  G(p^,u),t)  xn  Appendix  I. 

We  now  turn  our  attention  to  T.  ■*  ■+  _  ■*  +  ,  the  electron-atom 

kfrPffP*  l'Pi'Pl 

scattering  matrix  element.  This  is  a  well-studied  problem  which  we 
will  not  pursue  in  detail.16  Under  ordinary  situations  where  only 
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P  >>  Pf  leads  to  significant  cross-sections  in  Eq.(2.1),  the  Born 
approximation  can  usually  be  applied.  In  this  case 


P  -f  *►  ■+•  f  3  *•  CT  •  T*  i  i 

Ef ,Pf ,p;Ei,pi,p1  =  j  d  re'q  <f|v(r-x)|i> 


(2.30) 


where  fiq,  the  electronic  momentum  transfer,  is  given  by 


ftq  3  Px-  (m/mg)^  -  -  (m/mB)pf]  (2.31) 

— ^  ^ 

and  V(r-x)  is  the  electrostatic  interaction  between  an  electron 
(with  coordinate  r)  and  the  electrons  in  atom  B  (with  collective 
coordinates  x) . 


III.  Development  of  the  Green's  Function  in  the  Presence  of  an 
External  Field  -  Time-dependent  Spectrum 
We  wish  to  calculate  the  space-time  Green's  function 

G (x ' , t ' ; x , t) a  -i<T0|  T(aH(x'  ,t')ag(x,t)  1 1 4'0>  (3.1) 

(and  its  Fourier  Transform)  where  T  denotes  the  time-ordered  product, 
aH(x,t)  is  the  time-dependent  Heisenberg  picture  fermion  creation 
field  operator  and  | H,q>  satisfies  [cf.  Eq.(2.15)] 

<H(j  |4'0>  =  1-  (3*2) 

au(i,t)  is  the  Fourier  Transform  of  a„(p,t)[cf.  Eq.(2.14)]: 

n  n 

i±  + 

aH(x,t)  =  - — j  d3p  e  aH($,t).  (3.3) 

We  assume  that  the  separate  problems  of  the  adatom  orbitals 
<$k(x)  and  the  self-consistent  eigenstates  <}>CT(x)  of  the  unperturbed 
semi-infinite  surface  are  solved  so  that  the  time-independent  field 
operators  can  be  expressed  as 
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a(x)  =  [  ai4>i(x)  +  l  a0  <f>0  (x)  =  aR(x,0).  (3.4) 

The  surface  is  taken  to  be  metallic  and  a  includes  the  wave  number 

and  the  band  index,  a.  and  a  are  the  electron  annihilation 

1  a 

tu  i.u 

operators  for  the  i  n  adatom  state  and  a  surface  state,  respec¬ 
tively.  Each  electron  in  the  adatom-surface  plus  field  system  is 
considered  to  be  under  the  influence  of  a  self-consistent  field 
so  that  the  single-particle  Hamiltonian  can  be  written  as 

hAS(x,t)  =  h(x)  -  eEQ*x0  (t)  cosai^t.  (3.5) 

where  the  second  term  is  the  classical  interaction  Hamiltonian 
between  an  electron  at  x  and  an  external  laser  field  of  field 
strength  Eq  and  frequency  ioL  (in  the  dipole  approximation)  .  The 
field-free  Hamiltonian  h(x)  is  assumed  to  have  the  properties 

<(J)i|h]  <)>j>  =  (3.6) 

«t>alh|V>  =  ea5ao'  (3’7) 

where  and  e0  are  energies  belonging  to  the  states  <f>^  and  <|>a 

with  the  adatom  at  infinite  distance  from  the  surface.  The  interac¬ 
tion  leading  to  adsorption  is  given  by  the  matrix  elements 

< 4»i I h |  (t>0>  =  vi(y.  (3.8) 

Using  the  basis  {<p^,4>a},  h  can  then  be  written  in  matrix  form  as 


with  a,b  standing  for  {i}  or  {a}.  In  Eq.(3.11)  the  first  two 
terms  describe  "unperturbed"  electrons  corresponding  to  the  states 
|<J»i>  and  |<f>a>;  the  third  term  accounts  for  the  combined  effects  of 
adsorptive  and  radiative  interaction  between  the  states  |  >  and 

|  <J>  > ;  and  the  last  two  terms  describe  radiative  interactions  within 
the  sets  of  states  | $a>  and  respectively. 

Since  the  interaction  Hamiltonian  hj  only  involves  single¬ 
particle  coordinates  [cf.  Eq.(3.10)],  the  series  expansion  of  the 

Green's  function  [defined  in  Eq.(3.1)]  ssumes  a  particularly  simple 

21 

form.  The  usual  expansion  procedure  using  Wick’s  theorem  leads 
to  the  following  diagrammatic  representation  for  G{x',x): 


x' 


X' 

x' 

X  ' 

x’ 

1 

1 

»  1 

'  u 

t 

x,  1  ' 

i 

V  x2 

3 

A 

>  =  ' 

* 

+ 

X 

+  / 

,  +  X2 

t 

l  X.  1 

1 

4 

X1  " 

> 

1 

X 

1 

X 

>  < 

X 

1  X1 

(3.14) 


where  the  contracted  notation  x  stands  for  the  space-time  point 
(x,t) ,  etc.;  the  double  bar  and  the  single  bar  stand  for  the  full 
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Green's  function  and  the  "non-interacting"  Green’s  function 
respectively;  and  the  wavy  lines  stand  for  the  interactions 
hj  integrated  over  intermediate  space-time  variables.  Eq.(3.14) 
is  equivalent  to 

G(x'x)  =  G°(x',x)  +  ^  jd4xi  G°(x' ,x1)hI(x1)G0(x1,x) 

+  (^)  2  jd4Xid4x2G°  (x'  »x2)hI  (x2)G°  (x2  ,x1)hI  (x1)G°  (x^x) 

+  (%)  3|d4x1d4x2d4x3G°  (x'  ,x3)hI(x3)G°(x3,x2)hI(x2) 

X  G°(x2,x1)hI(x1)G°(x1,x)  + .  ,  (3.15) 

where  G°  is  the  "non-interacting"  Green’s  function;  and  it  leads 
to  the  Dyson's  equation 

G  (x '  ,  x)  -  G°(x',x)  +  Jd^G^x' ^Ufx-^GU^x)  (3.16) 

where  Z(x),  the  proper  self-energy,  is  simply  given  as 
h  (x) 

I(x)  »  .  (3.17) 

Equation  (3.16)  can  also  be  represented  in  the  diagrammatic  form 


V  '  X  * 
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where  the  shaded  circle  stands  for  the  proper  self-energy . 

We  are  now  in  a  position  to  calculate  G (ic '  ,ic;u>,t)  ,  the  restricted 
Fourier  Transform  of  G(x' ,t ' ;x,t) ,  by  converting  Eg. (3.16)  into  an 
algebraic  equation.  First  we  note  that  the  QSA  allows  us  to  replace 
Z(x1,t1)  in  Eq.(3.16)  by  £(x^,t)  [cf.  discussion  following  Eq.(2.28)], 
Dyson's  equation  can  then  be  written 

G(x',t';x,t)  =  G° (x' ,t‘ ;x,t)  +  jd^jJdtjG0  (x'  ,t'  jx^t^  E  (x^,t) 


X  Gfx^t^Xjt) . 

Fourier  transforming  with  respect  to  t'-t,  Eq. (3.19)  becomes 


(3.19) 


tG(x'  -  G°(x',x,id)] 


X  G  (x.  ,x  a) '  ,  t)  , 

F 


(3.20) 


since  G°  only  depends  on  the  difference  of  the  time  variables,  hQ(x) 
being  time-independent.  By  first  performing  the  integration  with 
respect  to  tj  and  then  w',  the  right  side  of  Eq.(3.20)  can  be  ex¬ 
pressed  as 

— i-jjdu)e~i(i)t  ’  Jd' 3x1G°  ( x ' ,  Xj_ ,  w)  l  ( xx ,  t)  fdw  ’ ela) '  fcG  ( , x , u '  ,  t) 

X  |dt1ei(M-“,)tl 

=»  _ijdu)e"iu)t '  |d3x1G°  (x'  ,xx  ,u>)  I  (Xj. ,  t)  e1U)tG  (^  ,x,u>  ,t) 

=  -i|dwe“iU>  (t_t' J  jd^j^G0  (x' ,xlfw)  E  (x1,t)G(x1,x,w,t)  .  (3.21) 


VV 
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Eq.(3.20)  then  leads  to  the  following  Dyson's  equation  for  G(x* ,x,w,t) : 

G(x',x,u>,t)  =  G°(x',x,ui)  +  jd3x^G^  (x‘  ,x^,<u) E  (x^, t)  G  (x^ ,x,w, t)  .  (3.22) 

Next,  the  spatial  Fourier  transforms  of  G(x',x,w,t)  are  performed 
to  obtain  G(icf  ,Jc,ai,t)  .  Eq.(3.22)  can  be  written  as 


-—«•  d3k'd3kelk' *x'  elk ‘x [G(jt '  ,£,u),t)  -  G°  (ic •  ,ic,w)  ] 


-i^|d3x1  d3k'd3k1e1^' ‘x'e1^1'xlZ(i<,,ic1,tl),t)|d3k2d3kelk2  ** 


x  eik  XG(lc2,ic,u>,t)  • 

where  Z  (ic '  ,oj , t)  =  jd3x'd3x^' 


^•:iki,xuftf  * 


G  (x'  ,x^,w)  E  (x^,t) 


(3.23) 


(3.24) 


Again,  by  doing  the  appropriate  integrations,  we  can  reduce  the 
right  side  of  Eq.(3.23)  as  follows: 

— d3k  'd3kelk'x ' e1^ '  x| d3kxZ  (£  ,1^ , u ,  t)  jd3k2G  (1<2  ,t  ,uj  ,  t) 


fd3x.e1  (kl+k2)  ,X1  =  — L_.fd3k’ 
>  1  (2ir)yJ 


.3.  ik’>x'  ik*x 
d  ke  e 


X  |d3k^Z0c' ,£^,u>,t)G(-k^,£,w,t)  . 


(3.25) 


Comparison  with  Eq.(3.23)  yields  the  Dyson's  equation  for  G(k',k,u,t): 


G(Jc' ,&,u),t)  =  G°(ic',ic,u)  +  — ^—y[d\.  Z  ()c  '  ,-jc.  ,a>,  t)  G  (Jc.  ,jc,a>,  t)  ,  (3.26) 

(2tt)  1  1 


which  can  also  be  written  in  the  diagrammatic  form 
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k' 


ft 


it' 


k- 


/f 


k* 


+ 


k 


(3.27) 


where  the  angular  bar  stands  for  the  interaction  Z/(2ir)^  and  the 
intermediate  momentum  variables  are  to  be  integrated  over. 

At  this  point  it  is  convenient  to  expand  G  and  G  as 


G (k  '  ,1c ,u ,  t)  =  l  <^*(-k’)<(>B(ic)GaB(aj,t), 

G°  ()c  1 ,  1c ,oj)  =  l  V-£’)<J>6  (£><£<<* >)5ag 
aB 

where  each  of  the  sums  runs  over  the  complete  set  of  indices  {i} 
and  {a}  [i.e.,  the  non-interacting  adatom  and  surface  states,  cf. 
Eq.  (3.4)J  and 


(3.28) 

(3.29) 


>„<*>  - \ 


d3xe-iS-;0„(J). 


(3.30) 


Expressing  G° (x ’ , x^ ,u)  also  as 

G°(x',x1,oi)  =  l  4>*(x')(|>ct(x1)G°  (u>)  , 


(3.31) 


we  obtain  the  following  form  for  Z(k’,k,u),t) 


’’SE 
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Z<i?\iV«.,t)  -hi  fd3x'd3x ($,) 

a  ■*  a  a  j. 

X  hjtXjrtjG^w)  =  l  ^(-^'JXo^ftJG^w)  ,  (3.32) 


where 


.  r-  -*• 
m  ^XJC  *  X 

Xa(ic/t)  =  i-jd3x1e  1  ■L4)a  ,t)  . 

Using  Eqs. (3.28,  (3.29),  (3.32)  and  (3.33),  Eq. (3.26)  becomes 

l  V-iS->i6<)E>tGo6<.>,t)  -  5°(»)6ae) 


a, 6 
1 


(3.33) 


(2ir) 

_ 1_ 

(  2it  ) 


a'p 


J,  fdVa<W<*i'Wd3kie  1  1 


*  ^a'(“^l,Ga'8(a,'t) 

- 1  ♦j(-i?',;6(i?)5“(u,)[i,jo,(t>aa,6(U,t). 


a8 
where 


Ya'a(t)  5  {d3x1*o'{$i)'S(*l't)*o(*1) 


(3.34) 


(3.35) 


may  be  considered  to  be  a  generalized  time-dependent  Rabi  frequency. 
Eq.(3.34)  yields  immediately  the  algebraic  Dyson's  equation  for 


(3.36) 
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or  J'1 =s“e‘ 

(3.37) 

In  matrix  notation, 

G(<o,t)  =  r_1(to,t). 

(3.38) 

where 

p  -  _  y*  (t) 

“6 

(3.39) 

Since 

G°(<o)  =  — —  , 
a  to-toa  ' 

(3.40) 

where 

“a  E  V*' 

(3.41) 

raB-  "V- 

(3.42) 

If  we  let 

"oB  m  V  ' 

(3.43) 

we  have 

V«  *  "S  V“ije(t)  oos\.t 

(3.44) 

*  S<via  '  ^0 ^ oosw^t) , 

(3.45) 

and 

W(t)  *  " •S^o'"<jo,9<t)  "n,1- 

(3.46) 

Eqs . (3.28)  and  (3.38)  then  determine  the  time-defendent  Green's 
function  G(ic'  ,& ,ui,t)  . 

G  is  related  to  G^,  the  advanced  Green's  function,  through  the 
following  relationships: 

G A ( p ,  to ,  t )  =  G  (p  ,p  ,io ,  t )  ,  u>  <  y(t)/ft 


=  G*(p,p,oj,t) 


oj  >  U(t)/fi 


(3.47) 
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where  y(t)  is  the  time-dependent  chemical  potential  of  the  adatom- 
surface-field  system.  According  to  the  Lehmann  representation, 22 
u(t)  is  determined  as  the  point  in  real  w  space  at  which  the  imagi¬ 
nary  parts  of  the  poles  of  G  (regarded  as  a  function  of  oj)  change 
sign.  These  poles  are  determined  by  Eqs.(3.38)  and  (3.42).  Eq.(2.26) 

can  then  be  used  to  give  the  spectral  function. 

IV.  Conclusion 

We  have  proposed  a  non-perturbative  formalism  to  treat  the  bound- 
continuum  (BC)  problem  of  collisional  ionization  of  solid  surfaces. 
This  formalism  has  certain  advantages  over  conventional  treatments 
of  BC  problems.  First  it  provides  a  direct  way  of  obtaining  differ¬ 
ential  cross  sections  with  respect  to  the  continuum  of  emitted  elec¬ 
tronic  energies,  whereas  the  usual  optical  potential  methods  would 
only  lead  to  the  total  cross  section.  Second,  it  bypasses  the 
approximation  of  limiting  the  problem  to  a  finite  number  of  chan¬ 
nels,  usually  two,  that  is  often  required  in  coupled-channels  ap¬ 
proaches.  Further,  it  renders  the  discretization  of  electronic 
energies  unnecessary,  which  is  again  the  standard  procedure  in 
coupled-channels  as  well  as  semiclassical  treatments. 

The  QSA  and  the  IMA,  however,  pose  limitations  on  the  validity 
and  usefulness  of  the  present  approach.  In  situations  where  the 
applied  field  has  frequencies  large  compared  to  or  of  the  order 
of  1/t ,  the  QSA  breaks  down,  although  the  IMA  may  still  be  applicable 
to  the  field-free  problem.  Fields  with  large  may  by  themselves 
lead  to  electronic  excitation  or  ionization,  and  the  inclusion  of 
collisional  effects  may  not  necessarily  provide  extra  information 
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on  the  systems  of  interest.  However,  the  general  collisional  problem 
with  arbitrarily  variable  w  is  definitely  of  immense  theoretical 
interest.  It  should  also  be  noted  that  the  underlying  criterion  for 
the  validity  of  both  the  QSA  and  the  IMA  is  the  shortness  of  the 
collision  time  t.  With  the  relaxation  of  this  criterion  to  the 
extent  that  the  IMA  is  no  longer  valid,  electron  correlation  effects 
will  have  to  be  considered  which  would  render  the  Green's  function 
formalism  much  more  complicated  than  the  present  treatment. 
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In  this  appendix  we  will  derive  the  Dyson's  equation  for 
the  quantity  G(p^,u),t)  [cf.  Eq.(2.29)].  We  start  f  rom  Eq .  ( 3 .16)  : 

G(x',t';x,t)  =  G°  (x' ,t '  ;x,t)  +  d3x1  dt^0  (x '  ,  t '  ;  x^  ,  ^ ) 


X  Z  (x^t^Gfx^t^-Xjt) 


(A  •  1 ) 


It  should  be  noted  this  equation  differs  from  Eq.(3.19)  [which  is 
used  to  obtain  the  restricted  Fourier  Transform  G(p^,w,t)]  only 
in  the  time  variable  in  the  proper  self-energy  Z.  Fourier  Trans¬ 
forming  with  respect  to  t'-t,  Eq.(A.l)  becomes 


■ji-  dwe  [G  (x' ,x,w,t)  -  G°(x',x,u))] 


— — j  jd3x^  dt^Z  (x^  ,  t^)  | due 


-io)(t'-t,)  n  r  -iw'(t,-t) 

Gu (x' ,x^ ,w) I dw'e 


X  G(x1,x,oj' ,t) 
From  Eqs.(3.17)  and  (3.10) 


(A. 2) 


ZCx^t^  =  ^{h^txj^)  -  eE0  •x16  (t1)  cosajj^tj^} . 

The  t^-integration  on  the  right  side  of  Eq.(A.2)  can  then  be 


written  as: 


dt1Z (x1 ,t^) e 


-i  (a)  ’  —oj )  t. 


=  ^  h1(x1)  6  (w'-oj)  -  ieEQ-j^  j  dtj.  e'1  (tiJ ''cu)  fcl0  (tL)  cosa>Lt1 


h  h1(x1)6  (w'-w)  -  -^eEg  xi^(jJ+UL-a)'+in  +  u)-u>L-a)k+in}  , 


(A. 3) 
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where  n  is  a  vanishingly  small  positive  quantity,  and  the  inte¬ 
gral  representation  for  0(t)  [Eq.(2.23)]  has  been  used.  Perform¬ 

ing  the  (o'  integration  next,  we  find  that 

G(x',x,u»,t)  -  G°(x,,x,w) 


f  3  -  -  h.tx.)  ^  _  et;0‘xl 

j  d  x^Gg (x ' , x^,o))[— — — — G  (x^  , x ,o) ,  t)  -  — ^ — 

iu)  tf,  „  iu)"t 

x  {e  j§¥r  G(xi'x'a,"+(J+wL't) 

-ia)_  t  .  „  ia)"t 

+  e  2ttT  ai " -  in  G  (Xx  'x/ w”+a)-u>L, t)  }  ] 


eEA-x. 


-ia)Lt  I 


(A. 4) 


The  spatial  Fourier  Transforms  are  then  performed  by  rewritincr 


Eq . (A . 4 )  as  follows: 

1  fj3.  '3.  ik'-x'  ik  •  x  ,  ~  ,  *■  ..  ~0  ,,  ,  .  .. 

- d  k  d  ke  e  {G (k ' ,k , w, t)  -  G  (k  ,k,w)j 

(  2tt  )  6  J 


|d3x1  d3k'd3k1e1^' ’x'e  1  1Z1(k,,k1,w)  d3k2d3ke  2  1elk'x 


id)  t 

w  t  .  1  e  u  f  ,3  f  ^3,  ,,3,  ik-x'  iKl  X1 

X  C(k_,k,w,t)  +  - r-y  — ~ —  d  x..  d  k'd  k.e  e 

(2v) 

'“’fir  5^fdV3kelk2 


.  -*■  ik  •  x. 

ik-x’  1K1  X1 
i  e  e 


id)  "t  i 


1  e"1U)LtfJ 3  [^3,  , j 3,  iic'-x’  lkl’xl„ 

t^^-)d,!l)dkdkie  e  ^ 


t  xfdd)" 

,kl,u,)  J  2*1 


•  ,  + 

X  %,  t  ■  fd3k,d3ke  2  1elk  ’ XG  (k- ,k  ,d)"+d)-d)  ,t) 
0)  -in)  2  2  L 


where 


T  3  3  -j  J  i  ,y '  — iki  ‘ x.  h.  o  ->  -*• 

Z1(k’  ,^,00)  =  d3x'd3x1e  lk  x  e  1  —  (x^  Gu  (x '  ,Xj_,  w)  , 


(A. 5) 


(A. 6) 
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Z2(£«,£lfa))  =  |d3x'd3x1i 


-i^’-x'  "ikl'xl  eE0*xl  o 

IK  x  1  i- - u — -U^U  (Z< 


£ - i)Gu  (x' /Xl,w)  .  (A. 7) 


By  first  performing  the  x^  and  then  the  k2  integration,  we  obtain  the 
Dyson's  equation  for  G  {&'  ,ic,u),t)  : 


G(k' ,ic,u>,t)  -  G°  (5'  ,£,*»)  =  — L^.fd3k,Z.  (ic*  ,-H.  ,u)G(Z.1i,ta,t) 


.  „  iw"t  ~ 

— j  d3k2Z2  (k'  ,-k^cu)  [  ^  j  2iT  w"-in  G  (kl ,k  t) 


,  e -lu)Lt|'dto "  e^  ^  ,?■  t  ^  , 

2  JIttT  a)"r,-in  ff(kl'k,U)  +o)-u)L,t)]. 


(A. 8) i 


We  again  invoke  the  expansions  [cf.  Eqs.(3.28),  (3.29)  and  (3.31)]: 


G°  (x'  , x1  ,o>)  =  l  <t>*(x')4)a{x1)G°  (w)  , 


(A. 9) 


G°(k\k,u))  =  ^*(-it')i6(^)G°(a))<5aB, 


(A.  10)1 


G(  k',Jc,w,t)  =  l  ♦  *(-£')♦*  (£)<?„«(«, t) 
a6  a  6  a6 


(A. 11) 


The  interaction  terms  Z^  and  Z2  can  then  be  written: 

+  _  f  3  3  —  i  Tr  ’  •  v  1  *  X*l 


->  rf  3  3  -i]c 1  •  x 1  1ki*xt  -4.  hi  "o 

Z1(k',-k1,u)  =  [jd  x'd3xLe  1K  x  e  1  x»*  (x  ’ )  <t>a  (x±)  ■■  G  ^(u) 


=  I  ♦a(_k')xa1)  ("k!)Ga(w)' 


(A. 12) 


Z2(k’,-k1,aJ)  =  l  fd3x'd3x1e’ik  ‘x  e  1  1<j>*  (x' )  <}>a  (x^^)  ( - 2fT~)  Ga  ((o) 

cr 


,'eV*i  -0 


(A. 13) 
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Figure  1.  Schematic  representation  of  the  transition  matrix 
element  T(t).  The  shaded  circle  represents  the 
hole  momentum  distribution  function  v  (+)  ...  i  2  ,-*•  x 

<  '  o  I  a 

for  the  adatom-surface  +  field  system  whereas 

the  blank  circle  represents  the  electron-atom 

scattering  matric  elementT-  -*■  ±  _  •+•  ■*  .  At 

Ef ,pf ,p;Ei,pi,p1 

the  shaded  circle,  energy  is  adiabatically  conserved? 
and  at  the  blank  circle,  energy  and  momentum  are 
strictly  conserved. 
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